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We discuss the composite fermion pairing state in bilayer quantum Hall systems. After the 
evaluation of the range of the inter-layer separation in which the quantum Hall state is stabilized, 
we discuss the effect of inter-layer tunneling on the composite fermion pairing state at v — 1/2. We 
show that there is a cusp at the transition point between the Halperin (3, 3, 1) state and the Pfaffian 
state. 
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I. INTRODUCTION 

As a remarkable example of strongly correlated elec- 
tron systems, the fractional quantum Hall effect has been 
studied extensively. In the fractional quantum Hall 
systems, the Coulomb interaction plays a dominant role. 
Due to the presence of the strong magnetic field, the 
Coulomb interaction gives rise to a two-body correlation 
of non-zero relative angular momentum. The Laughlin 
wave function, Q which captures the essential properties 
of the system, consists of this two-body correlation only. 
The role of the Coulomb interaction to generate the two- 
body correlation becomes clear when it is described in 
terms of Haldane's pseudopotential. Q From the analysis 
of the pseudopotentials, one can see that the most funda- 
mental contribution comes from the short-range Coulomb 
interaction. 

In the single-layer quantum Hall systems, it is hard to 
change the short-range Coulomb interaction. Whereas 
in the bilayer quantum Hall systems, we can adjust the 
inter-layer Coulomb interaction by varying the inter-layer 
separation d. In addition, there is another parameter, the 
inter-layer tunneling amplitude t to control the system. 
Since there are such additional parameters to control the 
system, we can expect rich phase diagram in the bilayer 
quantum Hall systems. 

In the absence of inter-layer tunneling, the properties 
of the system seems to be well described by the Halperin 
(m,m,n) wave function which is an extention of the 
Laughlin wave function: 
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where f and J, are the indices for layers and z° = x° — iyj 
is the complex coordinate of the j-th electron. The inte- 
ger numbers m and n is associated with the total Landau 
level filling v as v = 2/(m + n). 

At v = 1/2, the (3,3, 1) wave function has good over- 
lap with the wave function of the finite size system m in 



the absence of the inter-layer tunneling and in the appro- 
priate region of the inter-layer separation. It was pointed 
out by Ho that the (3, 3, 1) state and the Pfaffian state 
H , which may be stabilized in the strong inter-layer tun- 
neling limit, H are unified to the form of the p-wave 
pairing state of composite fermions. 

In this paper, we discuss the stability of the composite 
fermion paring states in the bilayer quantum Hall sys- 
tems. First, we discuss the range of the inter-layer sep- 
aration d in which the (to, to, n) state is stabilized and 
the appropriate choice of the number of fluxes for com- 
posite fermions. After deriving the mean field equations 
for the composite fermion pairing states, we consider the 
effect of the inter-layer tunneling t on the p-wave pairing 
state at v = 1/2. The evolution of the (3,3,1) state to 
the Pfaffian state is shown being based on the analysis 
of the gap equation. We also examine the possibility of 
Haldane-Rezayi state |10 at v = 1/2. 

The outline of this paper is as follows. In Sec. [nj, we de- 
termine the appropriate choice of (to, n) and we describe 
the relationship between the Halperin (to, to, n) state and 
the p-wave pairing state of composite fermions. The for- 
mulation for the analysis of the pairing state is given in 
Sec. III. The gap equations for the triplet pairing state 



and the singlet pairing state are derived in the presence 
of the inter-layer tunneling. In Sec. IV, we discuss the 
inter-layer tunneling effect. Finally, we summarize the 
results in Sec. |v|. 



II. TWO-BODY CORRELATION 

In this section, we discuss the bilayer quantum Hall 
systems in the absence of the inter-layer tunneling in or- 
der to find the appropriate choice of the number of at- 
tached fluxes for composite fermions. For the determina- 
tion of those numbers, the effect of inter-layer tunneling 
may be negligible because they are associated with the 
two-body correlation due to the short-range Coulomb in- 
teraction as discussed in Introduction. 

Since the two-body correlations, which is connected 
with the numbers to and n in the (to, to, n) wave func- 
tion (Fy), is associated with the short-range Coulomb 
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interaction, m and n are determined by "high-enrgy" 
physics. We compare Haldane's pseudopotentials with 
various choices of to and n |ll[] because "high-energy" 
physics is governed by Haldane's pseudopotentials. 

The basis for the two-body electron correlations is 
given by the wave function for an electron pair with the 
relative angular momentum to and the angular momen- 
tum of the central motion being zero: 
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where £b = \J ch/eB is the magnetic length. Since Hal- 
dane's pseudopotential is equivalent to the Coulomb en- 
ergy estimated in first order, the total energy for "high- 
energy" physics for total N electrons is given by 
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with e the dielectric constant and A = d/2£g. In the 
thermodynamic limit, we obtain 



E { c\m,n)/N' z = e(m, d = 0) + e(n, d). 
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Although the two-body correlation energies are mani- 
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festly overestimated in E c (m,n), the right hand side 
of Eq. (H) may be reliable as far as we are concerned 
with the pseudopotentials. For the choice of m and n, 
we cannot choose arbitrary pair of to and n. There is a 
constraint on the choice of to and n. From the Halperin 
(to, to, n) wave function, one can see that the angular 
momentum of the electron at the edge of the sample is 
equal to (N — 1) x m + N x n = M. Since the wave 
function of this electron is proportional to z M e 



the density of it has its maximum at r = \/2M£b = R- 
Of course irR 2 is the area of the system. Taking the 
thermodynamic limit N — > oo, we obtain 
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We consider the case of symmetric electron density pi = 
Pl- In this case, Eq. (||) is reduced to 

m + n = 2/v. (7) 

Now we determine to and n which gives the lowest 
Ec(m,n) under the constraint ([?]). For the case of 



v = 1/2, the constraint (|7|) is to + n = 4. There- 
fore, the possible choice of (m,n) is, (4,0), (3,1) and 
(2,2). The pair (m, n) with m < n always has larger 
energy than that with to > n. Note that in case of 
(m,n) = (even, even), the system is not the quantum 
Hall state. It is compressible liquid of composite fcrmions 
and the total wave function can not be determined by the 
two-body correlations only. 

In Fig. p], we show the energy {m,n) / N 2 for 
(m,n) = (4,0), (3,1), and (2,2). The region where 
the choice of (m,n) = (3,1) gives the lowest energy 
E (2) (m,n) is 0.789 < d/2£ B < 1.480. The Halperin 
(3, 3, 1) state is stabilized in this region. 

For the case of v = 1, the constraint (0) is m + n = 2. 
Therefore, the possible choice for (to, n) is (2, 0) and 
(1,1). In Fig.0, we show the energy E^ ^ (to, n)/N 2 
for (2,0) and (1,1). The region where the choice of 

(2) 

(m,n) = (1,1) gives the lowest energy E c [m,n) is 
< d/2£ B < 0.703(= d c ). The Halperin (1,1,1) state 
is stabilized in this region. Though the above estimation 
is crude, the critical value d c for v = 1 is close to 2£b 
that was obtained by Murphy et al. experimentally p2[ . 

(2) — 

In general, the estimation of E c (to, n) shows that 

(2) 

the pair (to, n) giving the lowest E)j'{m,n) is (2/v, 0) 
for d 3> Ib- As we decrease the value of d, it changes 
as (2/^,0) ->■ (2/f - 1,1) -> (1/u, l/v). In this 

sequence, the quantum Hall state is stable at (to, n) = 
(odd, odd), whereas the compressible state of composite 
fcrmions is stable at {m,n) — (even, even). 

Now we discuss the relationship between the (to, to, n) 
wave function and the p-wave pairing state of composite 
fermions. 0] As shown in Ref. JtJ, the wave function of 
the p-wave pairing state of composite fermions at v = 1/2 
is, in the second quantized form, 
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where is the creation operator of the electron at 

z with spin a. For the (3, 3, 1) state, \ 1S given by 

Note that in this case the p-wave pair- 
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ing wave function is the (1,1,-1) state. There- 
fore, in general the (m, to, n) state is described by the 

p-wave pairing state with \ - 
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of composite 



fermions with the number of attached fluxes (4>i, 4>2) — 
(to — l,n + 1). (The even integer <f>i (^2) is for the 
intra(inter)-layer correlations.) For the Pfaffian state, % 
l/y/2 l/y/2' 
l/y/2 l/y/2 _ 

quantum Hall state, then the wave function of composite 
fermions is not the form of the pairing state. 
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III. COMPOSITE FERMION PAIRING 

In the last section, we have discussed the appropritate 
choice of the number of attached fluxes for composite 
fcrmions. Now we introduce composite fermions in the 
second quantized form: 



^(r)=^(r)e^«, 
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where the function J a (r) is given by 

J a {v) = Y J K a0 I d 2 r'pp(r')Imln{z-z'). (10) 

13 J 

Here p a (r) = ^£(r)^ a (r) = $+(r)$ a (r) and K = 
( ± j! 2 1 with (f>i and 4>2 being even integer. In terms 
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of the composite fermion fields and if), the kinetic en- 
ergy term of the electrons is rewritten as 
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Here a Q is the Chern-Simons gauge field, 



we restrict ourselves to the range of the inter-layer sep- 
aration, d where the states based on composite fermions 
are stabilized, we may neglect the three-body interaction 
term. In addition, we neglect the long-range Coulomb 
interaction, which gives rise to a pair-breaking effect, be- 
cause the pairing state of compoisite fermions may be 
stable in the region where the Halperin (m, m, n) state 
is stable. In the following analysis, we concentrate on 
the analysis of the pairing interaction and the inter-layer 
tunneling. 

Including the inter-layer tunneling effect, H t — 
—t J d 2 r ^|(r)^(r) + ip'^(r)ip^(r) , the Hamiltonian for 
composite fermions may be written as 
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where £ k = £ k 



-t. In 



2 /2m - (j, and £ k = ^ = 
the interaction term, we have restricted ourselves to the 
scattering processes of pairs with zero total momentum. 
Note that the formulation in this section can be easily 
extended to the multicomponet systems. 

From Eq.(15), we can define the mean field Hamilto- 
nian as 
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The Chern-Simons gauge field a a obeys the constraint 

V x a fv (r) =4> Y,K a p Pp (r), (13) 





where (j>o = ch/\e\ is the flux quantum. 

The first order term of Eq. (|TT| ) with respect to the fluc- 
tuation of the Chern-Simons gauge field A+a Q yields the 
minimal coupling term. Eliminating the Chern-Simons 
gauge field fluctuations upon using the constraint (13), 
we obtain 
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ki-k2\ 2 • This interaction gives rise 



where V Ulk2 = *g 
to an attractive interaction that leads to the p-wave pair 
ing state. jl3) The second order term of Eq. (11) with 
respect to the fluctuations of the Chern-Simons gauge 
field yields the three-body interaction term after elimi- 
nating the Chern-Simons gauge field fluctuations. From 
the analysis of non-unitary transformation, this three- 
body interaction term turns out to be the counter term 
to the short-range Coulomb interaction. |l4| However, if 
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where the pairing matrix is defined as 
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First we consider the triplet pairing case. Since we 
consider the symmetric bilayer systems, we take the sym- 
metric form of the pairing matrix: A k | = A k | and 
A k ^- = A k ^. Diagonalization of the mean field Hamil- 
tonian yields the following gap equations at zero temper- 
ature: 
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where = J(£ k - t) 2 + |A k T ± A k | 2 . Meanwhile for 

A k , the gap 

0{Ey - 1) 



the singlet pairing state, A 
equation is given by 
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where E^ = + |Ak| 2 . At zero temperature, this 
equation is reduced to 
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Note that there is the constraint E\j > t in the summa- 
tion over k'-space. 

In the absence of the inter-layer tunneling, we can take 
because a pairing state with A k ^ ^ may be 
stable only in the vicinity of the sample boundary. By 
this choice of the pairing matrix, the gap equation has 
the same form both for the triplet pairing state and for 
the singlet pairing state: 
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We can solve this gap equation by taking the form of the 
gap as Ak = AeFf(k, Uf) exp(— i£9k), where f(k,kp) = 
(k/kp) e for k < kp and f(k,kp) = (kp/k) 1 for k > kp 
and £ is an integer. |p|,[l4| From the analysis of this gap 
equation, we find that the ground state is the p-wave 
pairing state. Note that the s-wave pairing state is ex- 
cluded because the pairing interaction originates from the 
Lorentz interaction induced by the Chern-Simons gauge 
field interaction. [l4|] At v = 1/2, the p-wave pairing 
state corresponds to the (3, 3, 1) state. At v — 1 with the 
choice of {4>i,4>2) = (0, 2), the p-wave state corresponds 
to the (1,1,1) state. |j 



IV. EFFECT OF INTER-LAYER TUNNELING 



9 = 7r/4, we have the Pfaffian state. In Fig. ||, we show 
the inter-layer tunneling dependence of 9. Note that the 
Pfaffian state is stabilized in the r > 2 region. 

In Fig. ||, we show the inter-layer tunneling dependence 
of the gap A. Note that there is a cusp at r = 2. Re- 
flecting the presence of the cusp in the gap, the ground 
state energy also has a cusp at r = 2. 

Now we discuss the possibility of the Haldane-Rezayi 
state, or d-wave pairing state. Since the s-wave pairing 
state is excluded as mentioned above, the next leading 
singlet pairing state is the d-wave pairing state. From 
the analysis of the gap equation, we find that for the d- 
wave pairing state A ~ 1.2, which is smaller than that 
for the p-wave pairing state, in the region of r < r c where 
t c ~ 0.8 and the d-wave pairing state is not stabilized in 
t > t c . However, in the above analysis we have neglected 
the effect of the long-range Coulomb interaction. Since 
the effect of it is expected to be larger for the p-wave 
pairing state than for the d-wave pairing state, it might 
be possible that the d-wave pairing state becomes stable 
due to the effect of the long-range Coulomb interaction. 
In addition, impurites may affect the p-wave pairing state 
more than the d-wave pairing state. 



V. CONCLUSION 

In this paper, we have discussed the region of the 
inter-layer separation where the p-wave pairing state is 
stabilized and the effect of the inter-layer tunneling at 
v = 1/2. The Pfaffian state is stable above the critical 
tunneling strength and there is a cusp at the transition 
point between the Pfaffian state and the state continu- 
ously connected with the (3, 3, 1) state. 
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Now let us take into account the inter-layer tunneling 
effect. We consider the effect of it on the p-wave pairing 
state at v = 1/2. The possibility of Haldane-Rezayi state 
is discussed later. As we have discussed in Sec. 0, the 
appropriate choice of (</>i, (f>%) is (2, 2). Since the number 
of attached fluxes <f>\ and 4>2 is symmetric, we may take 
A k | = Akd and A k ^ = Ak&. In order to discuss the 
evolution of the pairing state, we define an angle as 

9 = tan" 1 -. (23) 
b 

This angle 9 characterizes the pairing state. For the case 
of 9 = 0, we have the p-wave pairing state that corre- 
sponds to the (3, 3, 1) state. Meanwhile, for the case of 
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FIG. 1. The two-body correlation energy E^\m,n)/N 2 
for (m,n) = (4,0), (3,1), and (2,2) in units of e 2 /e£ B . 
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FIG. 3. The inter-layer tunneling dependence of 0. 
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FIG. 4. The inter-layer tunneling dependence of the gap A. 
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